A consistent model for ttN transition distribution amplitudes 
and backward pion electroproduction 

J. P. Lansberg 1 , B. Pire 2 , K. Semenov-Tian-Shansky 2 ' 3,4 , L. Szymanowski 5 
1 IPNO, Universite Paris-Sud, CNRS/IN2P3, 91406 Or say, France 
2 CPhT, Ecole Polytechnique, CNRS, 91128 Palaiseau, France 
3 LPT, Universite Paris-Sud, CNRS, 91404 Orsay, France 
4 IFPA, departement AGO, Universite de Liege, 4000 Liege, Belgium 
5 National Center for Nuclear Research (NCBJ), Warsaw, Poland 

Abstract 

The extension of the concept of generalized parton distributions leads to the introduction 
of baryon to meson transition distribution amplitudes (TDAs), non-diagonal matrix elements 
of the nonlocal three quark operator between a nucleon and a meson state. We present a 
general framework for modelling nucleon to pion (ttN) TDAs. Our main tool is the spectral 
representation for irN TDAs in terms of quadruple distributions. We propose a factorized Ansatz 
for quadruple distributions with input from the soft-pion theorem for ttN TDAs. The spectral 
representation is complemented with a D-term like contribution from the nucleon exchange 
in the cross channel. We then study backward pion electroproduction in the QCD collinear 
factorization approach in which the non-perturbative part of the amplitude involves ttN TDAs. 
Within our two component model for ttN TDAs we update previous leading-twist estimates of 
the unpolarized cross section. Finally, we compute the transverse target single spin asymmetry 
as a function of skewness. We find it to be sizable in the valence region and sensitive to the 
phenomenological input of our ttN TDA model. 



PACS numbers: 13.60.-r, 13.60.Le, 14.20.Dh 



1. INTRODUCTION 

The familiar collinear factorization theorem |IJ |2| for exclusive electroproduction of 
pions off nucleons 

e(k) + N( Pl ) -> ( T *(g) + N( Pl )) + e(k') e(k') + tt(^) + N'(p 2 ), (1) 

valid in the generalized Bjorken limit (large Q 2 = —q 2 and s = (p± + q) 2 ; xb = and 
skewness variable £ = — being fixed 1 ; and small —t = (p2 — Pi) 2 ) gives rise to the 

description of this reaction in terms of the generalized parton distributions (GPDs) (see 
left panel of Fig. [D . 

According to a conjecture made in |3|, [|], a similar collinear factorization theorem for 
the reaction ([3]) should be valid in the following complementary kinematical regime: 

• large Q 2 and s; 

• fixed Xb and skewness variable £, which is now defined with respect to the w-channel 
momentum transfer: 

[Ptt + Pi) ■ n 

• the w-channel momentum transfer squared u = {p w — Pi) 2 (rather than t) is small 
compared to Q 2 and s. 

Under these assumptions, the amplitude of the reaction (|l|) factorizes as it is shown on the 
right panel of Fig.[l] This requires the introduction of supplementary non-perturbative 
objects in addition to GPDs - nucleon to pion transition distribution amplitudes (tcN 
TDAs). Technically, nN TDAs are defined through the ttN matrix element of the tri- 
local three quark operator on the light cone @-0: 

6^(Am, A 2 n, X 3 n) 

= E^^^fAmjlAm; X n]%^(X 2 n)[X 2 n; X n}^ (X 3 n)[X 3 n; X n] . (3) 

1 n is the conventional light-cone vector occurring in the Sudakov decomposition of the relevant mo- 
menta. 
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Here a, (3, 7 stand for quark flavor indices and p, r, x denote the Dirac spinor indices. 
Antisymmetrization stands over the color group indices 012,3. Gauge links may be omitted 
in the light-like gauge A ■ n = 0. 





FIG. 1: Left: Collinear factorization for hard production of pions in the conventional 
hard meson production kinematics. Right: Collinear factorization for hard production 
of pions off nucleons in the backward kinematics. 



The detailed account of this approach is presented in Refs. [|TO]-0. Apart from the 
description of hard exclusive pion electroproduction off a nucleon in the backward region, 
the same non-perturbative objects appear in the collinear factorized description of differ- 
ent exclusive reactions. Prominent examples are baryon-antibaryon annihilation into a 
pion and a lepton pair in the forward and backward directions |T3|-|H)|]. 

The physical picture encoded in baryon to meson TDAs is conceptually close to that 
contained in baryon GPDs and baryon distribution amplitudes (DAs). Baryon to meson 
TDAs are matrix elements of a three quark operator (i.e. with baryonic number one) and 
characterize partonic correlations inside a baryon. This gives access to the momentum 
distribution of the baryonic number inside a nucleon. The same operator also defines 
the nucleon DA which can be seen as a limiting case of baryon to meson TDAs with the 
meson state replaced by the vacuum. In the language of the Fock state decomposition, 
baryon to meson TDAs are not restricted to the lowest Fock state as DAs. They rather 
probe the non-minimal Fock components with additional quark- ant iquark pair: 

(Nucleon) = |***) + |^*^; ^) + .... ; 

(Meson) = |M) + |M; M) + (4) 



For baryon to meson TDAs one may distinguish the Efremov-Radyushkin-Brodsky- 
Lepage (ERBL)-like domain in which all three momentum fractions of quarks are posi- 
tive and two kinds of Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP)-like regions in 
which either one or two momentum fractions of quarks are negative. On Fig. |2|, we show 
the interpretation of ttN TDAs in the ERBL-like and in DGLAP-1,11 region within the 
light-cone quark model [l6j . As one can see from Fig.^| (a) the ERBL part is probing the 
non-minimal Fock components of the nucleon wave function. In the DGLAP-II Fig. ^| (c) 
region one rather probes the non-minimal Fock components of the meson state, while in 
the DGLAP-I Fig. |2| (b) region there is a non-vanishing contribution of the minimal Fock 
states of baryon and meson. This interpretation, obviously, is justified only at a very low 
normalization scale. The evolution effects may significantly change it at higher scales. 

Xi x 2 x 3 



Xx x 2 x 3 



(a) 



(b) 



Xl x 2 x 3 




FIG. 2: Interpretation of nN TDAs within the light-cone quark model fll6fl . Small 
vertical arrows show the flow of the momentum, (a): Contribution in the ERBL region 
(all Xi are positive); (b): Contribution in the DGLAP I region (one of Xi is negative), 
(c): Contribution in the DGLAP II region (two Xi are negative). 



Similarly to GPDs |T7|^iT?| , by Fourier transforming ttN TDAs to the impact parameter 
space (A T — > b T ), one obtains additional insight on the nucleon structure in the transverse 
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plane. This allows one to perform the femto-photography of hadrons [^D] from a new 
perspective. In particular, there are hints that nN TDAs may be used as a tool to 
perform spatial imaging of the structure of the meson cloud of the nucleon. This point, 
which still awaits a detailed exploration, opens a fascinating window for the investigation 
of the various facets of the nucleon interior. 

Our paper is organized as follows: In Sec. |2|, we provide a short summary of the basic 
properties of irN TDAs. In Sec. |3], we review the spectral representation for nN TDAs 
and propose a factorized Ansatz for the quadruple distributions constrained by the soft 
pion theorem for ttN TDAs. We complement the spectral representation with a D-term 
like contribution and build a two component model for ttN TDAs. Sec. |] contains the 
details of calculation of 7*iV — > ttN amplitude in the backward regime. In Sec. |5|, we 
compute the unpolarized cross section and transverse target single spin asymmetry of 
backward tt + and tt° electroproduction off protons using our two component model for 
7rN TDAs. Many technical details are relegated to Appendices A-D. Our conclusions are 
presented in Sec. ffl. 



2. GENERAL PROPERTIES OF ttN TDAS 



At the leading twist-3, the parametrization of the Fourier transform of ttN matrix 
elements of the three-local light-cone quark operator (Q) can be written as 

e <E ^** A * (p - n Va(p,r)|6^(Ain, A 2 n, X s n)\N L ( Pl )) 
= S ( Xl+X2 + X3 _ 20 Y,(U)f^ 8 PT, x H^(x u x 2 ,x 3 ,Z,A 2 ;v 2 F ) , (5) 

s.f. 

where P = Pl + Pw is the average momentum and A = p n — pi is the w-channel momentum 
transfer. The spin-flavor (s.f.) sum in (|5|) stands over all independent flavor structures 
ifa)^ 1 and the Dirac structures s pTtX relevant at the leading twist; i (a) is the nuclon 
(pion) isotopic index. The invariant amplitudes which are often referred to as 

the leading twist ttN TDAs, are functions of the light-cone momentum fractions Xi (i = 
1,2,3), the skewness variable £ (0), the w-channel momentum-transfer squared A 2 and 
the factorization scale /j,f. 



4(P 
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For given isotopic contents (say proton to 7r° TDA), the parametrization ([5]) of twist-3 
irN TDAs involves 8 invariant functions V±\ , A^\ 4 ( see ^q- ( |A 1| ) ) . However, not 



all of them are independent. Taking into account the isotopic and permutation symmetries 
(see [0), one may check that in order to describe all isotopic channels of the reaction 
(0), it suffices to introduce eight independent ttN TDAs: four in the isospin-| channel 
and four in the isospin-| channel. This result is analogous to the case of leading twist 
nucleon DAs: initially, the parametrization || involves three proton DAs V p , A p and T p . 
However, due to the isotopic and permutation symmetries, these three functions may be 
expressed through the unique leading twist nucleon DA = V p — A p . Neutron DAs are 
expressed as {V n , A n , T n } = {-V p , -A p , —T p }. 

The support domain of baryon to meson TDAs in the light-cone momentum fractions x^ 
(5^. Xi = 2£) was established in It is given by the intersection of three stripes — 1+£ < 
Xi < 1 + £. Instead of dealing with three dependent light-cone momentum fractions Xj, 
one can switch to the independent variables. A convenient choice of independent variables 
is the use of the so-called quark-diquark coordinates |H| . Due to the symmetry of the 
support of baryon to meson TDAs under rotations by the ^ angle, there exist three 
equivalent choices of quark-diquark coordinates (i = {1, 2, 3}): 



1 3 

Wi =Xi-£; v j = - \] e ik i x k , (6) 



2 

k,i=i 

where em is the totally antisymmetric tensor. The support domain of baryon to meson 
TDAs in terms of quark-diquark coordinates can be parameterized as: 

- 1 < Wi < 1 ; - 1 + | C - £ | < vt < 1 - | C - & | , (7) 

where £ = 

As pointed out in Ulij , the scale dependence of ttN TDAs is described by the appropri- 
ate generalization of the ERBL/DGLAP evolution equations. Splitting functions in this 
case turn out to be much more complicated, as they include different pieces in different 
(ERBL-like and DGLAP-like) kinematical regions. 

Exactly as for the case of the usual parton distributions and GPDs, evolution of itN 
TDAs can also be treated in terms of renormalization of local operators corresponding 
to the Mellin moments of vrA^ TDAs in Xj. Evolution properties of the local operators 
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in question were extensively studied in connection with the scale dependence of nucleon 
DAs (see Refs. |22| . f23|). The conformal partial wave expansion of ttN TDAs over the 
conformal basis of the Appel polynomials or the Jacobi and Gegenbauer polynomials 
represents an alternative strategy for the parametrization of ttN TDAs in the spirit of the 
dual representation of GPDs [24, or complex conformal partial wave expansion [26|. 



Similarly to the GPD case, the polynomiality property in £ of the Mellin moments 
of tcN TDAs in the light-cone momentum fractions Xi is the direct consequence of the 
underlying Lorentz symmetry. For n\ + n 2 + n 3 = N the highest power of £ occurring 
in the (m, n 2 , n 3 )-th Mellin moment of Vff , Ajjjf, is iV + 1, while for T 3 ^f it is N. 
Consequently, TDAs V^, AJ 2 , T* 2 include an analogue of the D-term contribution [27], 
which generates the highest possible power of £. 

The most direct way to ensure both the polynomiality and the support properties for 
ttN TDAs is to employ the spectral representation in terms of quadruple distributions, 



generalizing the familiar Radyushkin's double distribution representation for GPDs [28 



3]J . In phenomenological applications, an inviting strategy, which proved to be successful 
in the case of GPDs, is to construct a factorized Ansatz for the corresponding spectral 
densities. However, contrary to the GPD case, 7rN TDAs lack a comprehensible forward 
limit, £ — » 0. This hampers the construction of the hypothetic factorized Ansatz for 
quadruple distributions with input at £ = as suggested in |lT| . 

In this paper, we build up a consistent model for tcN TDAs relying on their chiral 
properties. Chiral dynamics constrains ttN TDAs in the opposite limit, £ — > 1. Indeed, 
ttN TDAs are conceptually much related to pion-nucleon generalized distribution ampli- 
tudes (nN GDAs), which are defined through the cross-conjugated matrix element of the 
same three quark operator (|3|): 

(0|O^(Am, A 2 n, X 3 n)\N L ( Pl )n a (-p n )) . (8) 

A similar correspondence was previously established between pion GPDs and 2n GDAs 



27 , |32|| . For simplicity, let us consider the pion to be massless (m = 0). In this case the 
point £ = 1, A 2 = M 2 , where M stands for the nucleon mass, belongs simultaneously to 
both physical regions: that of tcN GDAs and that of irN TDAs (see discussion in ||12|| ). 



Moreover, it is at this very point that the soft-pion theorem [33]| applies for nN GDAs. 



As argued in p3|, RSI, this allows us to express nN GDAs at the threshold in terms of 



the nucleoli DAs V p , A p and T p . In the chiral limit, the soft-pion theorem for GDAs also 
constrains ttN TDAs similarly to the way []36] the soft-pion theorem |32| for 2n GDAs in 
the chiral limit links the isovector pion GPD at £ = 1, A 2 = to the pion DA ip n . 

In the chiral limit, the soft-pion theorem thus provides the normalization point for ttN 



TDAs. The explicit form of the soft-pion theorem for ttN TDAs was established in JT2 
In this paper, we use this information as input for realistic modelling of ttN TDAs based 
on the spectral representation in terms of quadruple distributions. 



3. SPECTRAL REPRESENTATION FOR irN TDAS, FACTORIZED ANSATZ 
FOR QUADRUPLE DISTRIBUTIONS AND D-TERM 

In this section, we propose a two component model for irN TDAs involving the following 
contributions: 

1. a spectral representation with input fixed at £ = 1 from the soft-pion theorem; 

2. a nucleon-pole exchange in the ^-channel which is a pure D-term like contribution 
complementary to the spectral representation. 



To do so, we formulate the spectral representation constructed in |TT| in a way suitable 
for the implementation of chiral constraints for ttN TDAs. This allows us to propose a 
factorized Ansatz for quadruple distributions with input from the soft-pion theorem. 

1. Toy exercise: GPD case 

To give an idea of the new type of factorized Ansatz for quadruple distributions, let 
us first consider how one can constrain a GPD model from the £ = 1 limit rather than 
from the forward limit £ = 0. Let us consider the standard double distribution (DD) 



representation for GPDs ]2q pi 



H(x,g)= [ d(3da5(x- (3 -a£)f(p, a), (9) 
Jn 

where Q is the usual domain in the spectral parameter space 

fi = {|/3|<l; \a\<l-\P\}, (10) 



and f([3,a) is the double distribution. Let us perform the change of variables: a = 
This gives: 



h(x, o = J dn j des^x + ^e-^^ f(k, e) , (ii) 

where F(k, 9) = fO^, ^2)' Instead of the usual factorized Ansatz in ([]) in the 
variables {/3, a}, 

f(P,a) = h(P,a)q{P), (12) 
let us employ in ( p]) the following factorized Ansatz in the variables {k, 9}: 

F(k, 9) = 2<p(K)h(6), (13) 
with the profile h(9) normalized according to 

d9h{9) = 1 . (14) 



-i 

Obviously, Eq. (]TTf) then gives H{x^ = 1) = <p(x). 

In order to implement the so-called "Munich symmetry" [|37) f((3,a) = f((3, —a), 
which is the consequence of hermiticity and time-reversal invariance of non- forward matrix 
element entering the definition of GPDs, one should require that 



h(0) = <p(- 



i i - 1 

d9ip(-9) 

i 



(15) 



Let us emphasize that, in the GPD case, symmetry requirements unambiguously fix the 
shape of the profile h in the factorized Ansatz (|l~3|). 

Although it is but a toy model, the factorized Ansatz (0) may be applied to the case 
of quark isovector GPD of a pion H%~ d , which in the soft-pion limit |32j reduces to the 
pion DA ip n : 

hm#r d (x,£,t = 0) = ^(x). (16) 
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2. An alternative form of the spectral representation for irN TDAs 



Let us now apply the trick of previous subsection for the case of irN TDAs. According 
to the spectral representation for tiN TDAs reads: 

H^ N \ Xl ,x 2 , x 3 = 2£-x 1 -x 2 , 



8(xi -£-pi- a 1 £)5(x 2 - f - 02 - a 2 C) 



x8(/3i + P 2 + 03)8(a 1 + a 2 + a 3 + l)/(/3i, 02, 03, "i, "2, a 3 ), 



(17) 



where Qi denote three copies of the usual domain ([10]) in the spectral parameter space. 
The spectral density / is an arbitrary function of six variables, which are subject to two 
constraints imposed by the two last delta functions in eq. (|TT|). Therefore, / is effectively a 
quadruple distribution. The spectral representation flT7D by the very construction ensures 
the polynomiality and the support properties of irN TDAs. 

Let us employ the particular choice of the quark-diquark coordinates (wi, Vi) @ and 
introduce the following combinations of the spectral parameters: 

1 3 

Ki = ai + 0i] 9j = - ^ e ik i(a k + 0k); 

k,l=l 
3 



Hi = aci - 0i ; Xj = ^ e m {a k - k ). 



kl=l 



The spectral representation (fl~7| ) can then be rewritten as: 



-i 



2 



(i -£)-«*£) 



x<> f >k ~ 1 Xl (1 - - J" F i(«<> ^ A*i> A 



(19) 



The working formulas for the calculation of 7riV TDAs in the ERBL-like and DGLAP-like 



domains are summarized in Appendix [B2 . 

We suggest using the following factorized Ansatz for the quadruple distribution F, in 



Fifa, 9i, n h Xi) = 0i) h(n h Xi), 



(20) 
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with the profile function h(fii, Aj) normalized as 

,1 - Hjj 



d/Xj / dAj /i(/i i5 Aj) = 1. (21) 
-i ■/-^ 

Note that the support of the profile function h is that of a baryon DA. 

With the help of the spectral representation ([19|) , one can check that in the limit £ — > 1 
iYj now reduces to 

ffK u<, £ = 1) = V( Wi , Vi ). (22) 

We also note that the support properties of Fj(/%, 0j, /ij, Aj) in the 0j)-plane corre- 
spond to that of baryon DAs. It is thus natural to use the combination of baryon DAs 
to which nN TDA reduces in the limit £ — >■ 1 due to the soft-pion theorem as input for 
V{w h v^. 

Let us denote the combination of nucleon DAs, to which nN TDA H reduces in the limit 
£ — > 1, as V(yi,y2,y3). It is the function of three momentum fractions yi (0 < yi < 1) 
satisfying the condition YliVi = Then, according to the particular choice of quark- 
diquark coordinates in fll7|), one has to employ in (|20|) : 



2 ' 


4 


4 


1 - «a - 20 2 


k 2 + 1 


1 - « 2 + 20 2 


4 


2 ' 


4 


1 - « 3 + 20 3 


1 - «3 


-20 3 «3 + l 


4 


4 


' 2 



1 + l l-K!+20! l-Ki-201 

V(k 2 ,9 2 ) = jV 

V(k 3 , 3 ) = V ^ - - ^ j • (23) 

The profile function h(fii, Aj) also has the support of a baryon DA: — 1 < /ij < 1; 
< Aj < ±=p. Contrary to the GPD case, no symmetry constraint from hermiticity 
and time-reversal invariance occurs for quadruple distributions. Therefore, we are free to 
employ an arbitrary shape for the profile function. For example, we may assume that it 
is determined by the asymptotic form of a baryon DA (120yiy 2 y3 with J2i Hi = 1) : 

h(jH, A,) = ^ (1 + /i 4 )((l - /i,) 2 - 4A 2 ) . (24) 

In fact, this is the simplest possible choice for the profile vanishing at the borders of its 



domain of definition. The normalization (21) is obviously ensured. 



On Fig. [3], we present the result of the calculation of ir°p TDAs from the factorized 
Ansatz ( PD| ) with the profile function ( p^ ) as a function of quark-diquark coordinates 
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w = W3, v = t>3 defined in ([]). In accordance with the soft-pion theorem, in the £ = 1 
limit, our ir p TDAs p , p and v are reduced to the following combination of 
nucleon DAs JT2| 2 : 

if'fo, x 2) s 3 , £ = 1) = \ x ^ (y , y , y ) • (25) 

We employ the Chernyak-Zhitnitsky (CZ) nucleon DA as the numerical input. 

Our spectral representation provides a lively Xi and £ dependence for ttN TDAs. How- 
ever, the suggestion of a reasonable A 2 dependence still remains an open question. The 
most straightforward solution would be, similarly to early attempts in the GPD case, to 
try a sort of factorized form of A 2 dependence for quadruple distributions (|2~0|): 



Fi(Ki, e u m, Xi, A 2 ) = 4V(K h Oi) h(ja, Ai) x G(A 2 ), (26) 

where G(A 2 ) is the ttN transition form factor of the local three quark operator 
0^(0,0,0) (|). This leads to a factorized A 2 -dependence for nN TDAs: 

H* N { Xi , £, A 2 ) = H" N ( Xi , £) x G(A 2 ). (27) 

The determination of G(A 2 ) goes beyond the scope of the present paper. Let us however 
note that such a factorized form of the A 2 dependence is known to be oversimplified and 
was much criticized in the GPD case (see e.g. discussion in fl38|). 



2 



Note that eq. (11) of 15 and eq. (19) of 10 contain a sign error for p as well as erroneous overall 



factors 2. This affects the numerical results of these papers. 
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FIG. 3: n p TDAs V£ p , A\ p and p computed in the model based on the factorized 
Ansatz fl20|) with the profile function ( pif) as functions of quark-diquark coordinates w = 
u> 3 , v = t>3. In the limit £ = 1, as required by the soft-pion theorem, TDAs are reduced to 



the appropriate combinations of nucleon DAs V' p , A p and T p (see eq. |25p. For £ < 1, 7rA r 
TDAs are obtained by "skewing" £ = 1 limit. CZ nucleon DAs ||T] are used as numerical 
input. 
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3. D-term-like nucleon pole contribution 



Similarly to the GPD case, ttN TDAs within the spectral representation ([H]) do not 
satisfy the polynomiality condition in its complete form. As it was pointed out in ||12|| , 



the spectral representation for ttN TDAs Vi^ N \ , has to be complemented by 

an analogue of the D-term. TDAs do not require adding the D-term. This D-term 

has a pure ERBL-like support in Xi and hence it contributes only to the real part of the 
backward pion electroproduction amplitude (pTCf). In this paper, we employ the simplest 
possible model for such a D-term which accounts for the contribution of the u-channel 
nucleon exchange into ttN TDAs computed in |12| . This model shares many common 
features with the pion pole model for the polarized nucleon GPD E suggested in Sec. 2.4 
of Ref. H (see Fig. |). 





FIG. 4: Left: pion pole exchange model for the polarised GPD E; lower and upper blobs 
depict pion DAs; the dashed circle contains a typical LO graph for the pion electromag- 
netic form factor in perturbative QCD; the rectangle contains the pion pole contribution 
into GPD. Right: nucleon pole exchange model for tcN TDAs; dashed circle contains a 
typical LO graph for the nucleon electromagnetic form factor in perturbative QCD; the 
rectangle contains the nucleon pole contribution into tcN TDAs. 



The explicit expression for the contribution of the w-channel nucleon exchange into the 
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isospin-i ttN TDAs was established in ||12|| : 

{V lf A 1 ,r 1 } (7rAr)l / 2 (x 1 ,x 2) x 3 ,^A : 



^(940) 



^(940) 



0. 



{V 2 , A 2 ,T 2 } (7riY)l/2 (x 1? x 2 ,X3,e,A 
where we employ the notation 

3 

©ERBL^l,^,^) = n^^ - Xk - 2< ^' 



= i{^ 1) A 1) T 1 } ( ^(x 1 , a ; 2 ,xs ) e ) A 2 ; 

Af(940) 2 1 ' 



(28) 



(29) 



fe=i 



M denotes the nucleon mass; f m is the pion weak decay constant and g n NN is the pion- 
nucleon phenomenological coupling (see e.g. pi| ). 

The nucleon pole contribution into 7c + p and 7r°p expressed through isospin-1/2 nN 



TDAs (H) reads: 



{V li2 , A h2 ,T h2 y + v 

N{940) 

= -V2{V lfl , A lj2 ,T lj2 y°*> 

{Ts A V +P = {Ts,i}*° P 

^(940) 



Af(940) 



^(940) 



V2{V 1>2 , A 1;2 ,T 1>2 }^ 



/2 



^(940) 



0. 



(30) 



4. CALCULATION OF 7 *iV -> vriV AMPLITUDE 

Within the factorized approach, the leading order (both in a s and 1/Q) amplitude for 
hard exclusive 7*iV — > ttN reaction in the backward region, .M* s , reads ||10|| : 

7 14 



M: 



,(47ra s )V 47r « em /Ir 1 



1S2 



Q 4 
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J J \ a= i a= s 



(31) 



-K S2 N 3 x d 3 yhJ2Ta + J2 T « 

J J \ a= i a= s 

where f n = 93 MeV is the pion weak decay constant and /jy ~ 5.2 • 10~ 3 GeV 2 is a 
constant, which determines the value of the dimensional nucleon wave function at the 
origin; a em ~ is the electromagnetic fine structure constant; and a s is the strong 
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coupling constant. The convolution integrals in Xi and in fl3~l"|) respectively stand over 
the supports of ttN TDAs and nucleon DAs in T a and T' a . The spin structures S and S' 
are defined as 

S^ S2 = U(p 2 , s 2 )e(X)j 5 U( Pl , si) ; < S2 s ^(pa, S2 )e(A)A r7 5 f/(pi, si), (32) 

where e(A) denotes the polarization vector of the virtual photon. We introduce the fol- 
lowing notations: 

/■!+£ rl+Z 

{Z,Z'}(£, A 2 ) = / dxi / rfx 2 / dx 3 + x 2 + x 3 - 20 
j-i+C ^-1+5 

/■l r 1 r 1 ( 7 u > 

x / d Vl \ dy 2 / dy 3 8{ Vl + y 2 + y 3 - 1) 2 ^{T a , T^} + ^{T Q , T^} 
Jo Jo \ a=1 a=8 



C = -i 



and rewrite (|3"TD as 



.(4vra s )V4vra em /^ 



54/, 



(33) 



A4 A = C— 



(34) 



The expressions for the coefficients T a and T' a for the 7*p — > 7r°p channel are presented 



in the Table I of Ref. [ lOfl - The result for *y*p — > 7t + n channel can be read from the same 
Table with the obvious changes: 

Qu ^ Qdi Qd ^ Qu i 

V p , A p , T p -» V n , A n , T n = -V p , -A p , -T p ; 

v l,2 > ^1,2 5 1 1,2,3,4 ~^ K l,2 ) ^1,2 5 J 1,2,3,4 • 



(35) 



We note that in Ref. [10 a somewhat inadequate parametrization of nN TDAs was em- 
ployed. Within this parametrization, ttN TDAs do not satisfy the polynomiality property 
in its simple form due to the appearance of kinematical singularities (see discussion in 
Ref. |r2|). In this paper, we adopt the parametrization suggested in IT3] in which poly- 
nomiality is explicit. The relation between the two parameterizations is given by eq. 



As we note, x« and ?/, dependencies in coefficients T a (T^) are factorized. One therefore 
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anticipates that the convolution integrals in eq. ( |31| ) have the following generic structure 3 : 

J d 3 xK a (xi,X2,Xz) combination of ttN TDAs(xi, X2, x 3 ) 

x J d 3 yR a (y 1 ,y 2 ,y 3 ) combination of N DAs(y 1 ,y 2 ,y 3 ) ■ (36) 

K a (x 1 , x 2 , X3) and R a (yi,y2,y3) refer to parts of the singular convolution kernel in T a 
(T^) depending on Xi and yi respectively. These can be read from the Table I of Ref. fT0|j . 

The convolution integrals in yi in ( |36| ) are similar to those occurring within the per- 
turbative description of the nucleon electromagnetic form factor. The convolutions with 
singular kernels R a (yi, yi, 2/3) do not generate any imaginary part since the nucleon DAs 
have purely ERBL support and vanish at the borders of their domain of definition. These 
integrals can be calculated in a straightforward way. 



On the contrary, the convolution integrals in Xi with ttN TDAs in (36) may, in principle, 
generate a nonzero imaginary part of the amplitude. nN TDAs, indeed, do not necessarily 
vanish on the cross over trajectories Xi = 0, separating ERBL-like and DGLAP-like 
regimes, as well as on the lines Xi = 2£. 

Switching to quark-diquark coordinates (|), one may show that the following types of 
convolution kernels K a occur in (|36|): 

1 1 



Kji(wi,Vi) 



Throughout the following discussion, we adopt the convention that the first sign in the 
(±, ±) index of a quantity corresponds to the one in the w ± £ denominator while the 
second sign corresponds to the one in the v ± £' denominator. 
Thus, we have to deal with only two types of integrals: 

/J ± ' ± )(£)=/ dw dv- ,H(w,v,£), (37) 



and 



Here a = 1, 14 should be understood as a label. No summation over repeating a is implied. 
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for which we have to develop a method of calculation. The integration in (|37| ) and (|3q) 
stands over the support ([?[) of nN TDAs in quark-diquark coordinates. 

Using the formulas summarized in Appendix we establish the expression for the real 
and imaginary parts of 1^' 



J -i \ w + J-i+\e,-(,'\ \ v± t,) 
Rel i r ±) (0=v[ 1 dw—^—vf 1 l? ^ dv—!-—H(w,v,OTK 2 H(Z,0,O; 

Im/j + ' ±) (0 = ^/ dw^—H(w,Te,0 + ^[ dv-t—H(-Z,v,t); 
J-i w + 4 y_! w±4 

Im/;-' ±) (0 = TvrP/ dw-!—H(w,Tt;',0-irv[ dv-H^v^). (39) 

Let us now consider the second type of integrals by rewriting it as: 

4"' ±} (0 = ^ (w _ l i + iQ)2 {t ™h( w , t?, o + jw (to, o I , (40) 

where we introduced the notation: 

J (± W)=W dv——H(w,v,Z). (41) 

Let us emphasize that in (f40|) we are dealing with convolution of the product of two 
generalized functions with the test function H(w ,v, £). In order to assign meaning to this 
ill defined expression as it is done in eq. (|I0"D, H(w,=f£',£) and J^(w,^) and their first 
derivatives in w should be continuous in the vicinity of w = £. One can check that these 
assumptions are justified by the use of the spectral representation ( |TTD with continuous 
input quadruple distributions vanishing at the borders of their domain of definition. 

We obtain the following contributions to the real and imaginary parts of the amplitude 
from d40|): 



Re/^'^O 



(w-0 (w-0 



J-i (w-0 V 



(42) 
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The formulas for the calculation of the real and the imaginary parts of //^(C) an d 
In'^iO m the model based on the factorized Ansatz for quadruple distributions (|20| ) 
with input from the soft-pion theorem at £ = 1 and with the profile function h given by 
are summarized in Appendix [D[ 

We are going now to present the results of calculation of X(£) in our composite model 
for n N TDAs of Sec. [| As it was already pointed out, the coefficients T a , T' a (|3"3"D, which 
can be read from the Table I of Ref. jlTJ, are defined with respect to the alternative 
parametrization of ttN TDAs. The relation between that parametrization and the one 
employed in the present paper is summarized in Appendix [A]. Within the parametrization 
of Ref. HlQ| , X(£) receives contributions only from ttN TDAs {Vx, Ax, Tx, ^4 } while X'(£) 



receives contributions only from nN TDAs {V 2 , A 2 , T 2 , T 3 }^. 

One may establish the following relations for the spectral part of 7tN TDA model based 
on the factorized Ansatz fl20|) with input from the soft-pion theorem: 



{Vi, A x , Tx}* N (xi, x 2 ,x 3 , £) | j^ ctralpart = {Vi,Ax,Txy«(xx,x 2 ,x 3 ,0\^SJZ- (43) 



and 



{V 2 ,A 2 ,T 2 ,T 3 ,T^ N (xx,x 2 ,x 3 ,0C^ pait = 0- (44) 
Now we consider the nucleon pole part fl28|) of the two component model for ttN TDAs 

{Vx, A, Ti r N (xx, x 2 , x 3 , 1 § 40) = ^{Vx, A, Tx} 7rN (x 1 , x 2 , x 3 , || paper ; 

{V 2 , A 2 y N { Xl , x 2 , x 3 , 1| 940) = {Vx, Ax, Ti r N (xx,x 2 , x 3 , |f ^ ^ 

(T 2 + T 3 r( Xl , x 2 , x 3 , 1 @ = 2T^( Xl ,x 2 , x 3 , 1 j|* paper . (45) 



Consequently, in our model, the following relation holds for the nucleon pole contribution 
into X and I': 

A2 )|n ( 940) = Tr| ReX '^ A2 )|N(940)- ( 46 ) 

On Fig. |5] we present the results in our model for the real and imaginary parts of X(£) 
for backward production of tt° (two upper panels) and tt + (two lower panels) off proton 
showing separately the spectral part, the pole part and their sum. The CZ phenomeno- 
logical solution || for the nucleon DA is used as the numerical input for our model. For 
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small £ (£ ^ 0.3 4-0.5), the real part is dominated by the contribution of the nucleon pole. 
The contribution of the spectral part to the real part becomes relatively more important 
for larger £. The nucleon pole contribution is purely real. The apparence of a significant 
imaginary part stemming from the spectral component of the model is a distinctive fea- 
ture of our approach. This is crucial for the non-vanishing of the transverse target single 
spin asymmetry for backward pion electroproduction discussed in Sec. [5| 



CZ input nucleon DA CZ input nucleon DA 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 5: Real and imaginary parts of X(£) for 7*p — > 7c°p and 7*p — > ir + n backward 
production as functions of £ computed in our two component model for tcN TDAs. Dashed 
lines: nucleon pole contribution into ReX(£); thin solid line: spectral representation with 
input from the soft-pion theorem; solid line: sum of two contributions (for the real part). 



5. UNPOLARIZED CROSS SECTION AND TRANSVERSE TARGET SINGLE 
SPIN ASYMMETRY FOR BACKWARD PION PRODUCTION 

Let us first specify our conventions for the backward pion electroproduction cross 
section. In the one photon exchange approximation, the unpolarized cross section of hard 
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leptoproduction of a pion off a nucleoli ([[]) can be decomposed as follows [[45| : 

dV _ a cm (s - M 2 ) 

dsdQ 2 d<pdt ~ 4(2vr) 2 (^) 2 M 2 Q 2 (l - e) 

x (— + e _ + e cos 2^— + cos ) , (47) 

where ip is the angle between the leptonic and hadronic planes; s = {p\ + q) 2 = W 2 and 
t = (p2 ~ Pi) 2 are the Mandelstam variables; fcg is the initial state electron energy in the 
laboratory (LAB) frame (beam energy), e is the polarization parameter of the virtual 
photon: 



1+2 w-*y + *y tan2 | 



(48) 



where k' is the energy of the final state electron in the LAB frame and 0% is the electron 
scattering angle in the LAB frame. 

Within the suggested factorization mechanism for backward pion leptoproduction, only 
the transverse cross section receives a contribution at the leading twist level. Using the 
explicit expression relating scattering amplitudes of leptoproduction to those for virtual 



photoproduction (eq. (2.12) of Ref. we express ^y 21 in the center of mass (CMS) 



system of the pion and final nucleon through 7*iV — > Nn helicity amplitudes M.* 
defined in (|31|): 

dE'dtle'dti^ 

= -I^S^E{>^ + W) + ...}^x(^ + ... 

V S1,S2 



(49) 



Here, Q e i is the differential solid angle for the scattered electron in the LAB frame; fl n 
is the differential solid angle of the produced pion in N'n CMS frame (see Fig. |6| for the 
definition of angular variables); by dots, we denote the subleading twist terms supressed 
by powers of 1/Q] A is the usual Mandelstam function 



A(x, y,z) — \J x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. (50) 



T is the virtual photon flux factor in Hand's convention |46] given by 



q em k' L s-M 2 1 

2tt 2 k% 2MQ 2 l-e' [ ' 
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FIG. 6: Kinematics of electroproduction of a pion off a nucleoli in the CMS frame of 7* 
nucleon. 



Our present goal is to establish the expression for the LO transverse cross section 
through the helicity amplitudes defined in We rewrite our formula (pip as 



M? 32 = C^U( P2 , 8 2 ) T H Ufa, si), (52) 

where 

r H = e(A) 75 X-e(A)^75X'. (53) 

Let us now square the amplitude and sum over the transverse polarizations of the virtual 
photon and over the spin of outgoing nucleon: 



\ M t\ 2 = \ C \ 2 ^ E Tr { fe + M^H^Y^iPi + M) 7o rU}- (54) 
Let us first consider the trace 

^Tr{(p 2 + M)r„(p 1 + M) 7o r5,7o} = (55) 

-^(\)^{X)Tr{fa + M){f^X - 7 ^ 75 X')(pi + M)( 75 7"X* - 75 ^(J')*)} . 

We employ the relation 

V e"(A)e^(A) = -g^ + T ^—(p fl n u + p u n») (56) 

VP-n) 
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to sum over the transverse polarizations of the virtual photon. We use the backward 
kinematics for the reaction ([[]) summarized in fllQfl : 



l + £ M 
Pi-n= pi -p 



2(1 + 



p 2 -n = 0(l/Q 2 ); P 2-P=^F + 0(Q°). (57) 



Then for the part which is independent of the nucleon spin we get: 



Tr{(p 2 + M)T H (p 1 + M) 7o r^ 7o 



A j 1 



= 2Q ! (l ± « |I|2 _2Q ! (l ±i )^ |I , |2 + o(Q0) (58) 

Now we turn to the nucleon spin dependent part of the trace (54). 
^Tr{(p 2 + M)Y Hlh hipi + M)7 r^7 } 

X'j' 

= ^e^(A)e^(A){(p 2 + M)Yh(Pi + M) 7 5^V}x(X')* 
+ £V (A)e"*(A){(p2 + Mh^hiPi + M) 75 7 M }X'(X)* 
= 4 g2( ^ e(n, p, Sl , A T ) ( - zX(X')* + iX'(X)*) 

= -4 g2(1 e +0l ^ l |^i| sin(^ - ^)Im(X'(X)*). (59) 

In the last line of fl59"D , we consider si as being purely transverse and choose the reference 
frame so that si has only an x-component. Then 

e(n,p, Si,A T ) = i|A T ||si| sin(v? - v? s ), (60) 

where </? is the angle between the leptonic and hadronic planes and </> s is the angle between 
the leptonic plane the transverse target spin (see Fig. ^[). We employ the conventions in 
which e 0123 = 1 with 75 = — ^£ >iupa 'y^'yu1p1a ■ Finally, we conclude that 

I-Mt 1 ! 2 = ( |X|2 " ^ |X? ~ 2 ^ l?l1 Im(x ' (xr) sin ^) + oWWV 

where (p = <p — ip s . 
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Hence, we establish the following formula for the LO unpolarized cross section 
through the coefficients X, X', introduced in (|34|): 



d 2 (JT = 2 1 A(s,m 2 ,M 2 ) l + e nx|2 _ ^T| Z /,2A ffi2l 
dfi* 1 1 Q 6 128vr 2 s(s - M 2 ) £ ^ M 2 ' 1 J 

Within our kinematics 

A 2 = ^ fT A^ UL , (63) 

where m is the pion mass. 

On Fig. ^ we present our estimates for the unpolarized cross section ^p: Q f 
backward production of 7r° and n + off protons for Q 2 = lOGeV 2 and u = — 0.5GeV 2 in 
nb/sr as function of Xb- 

We use the two component model for ttN TDAs presented in Sees. |3] and f|. In order 
to quantify the sensitivity of our model prediction on the input nucleon DAs, we show 
the cross section for the case of four phenomenological solutions fitting the nucleon elec- 
tromagnetic form factor: CZ ||] (solid lines), Chernyak-Ogloblin-Zhitnitsky (COZ) [[y| 
(dotted lines), King and Sachrajda (KS) |42| (dashed lines) and Gari and Stefanis (GS) 
|43j (dash-dotted lines). The magnitude of these cross sections is large enough for a de- 
tailed investigation to be carried at high luminosity experiments such as J-lab@12GeV 
and EIC. We recall that the scaling law for the cross section (pi) is 1/Q S - 
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CZ v.s. COZ v.s. KS v.s GS input nucleon DA 
lO.Or — ■ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — - 




N° 0.2 



0.1 



0.0 0.2 0.4 0.6 0.8 

CZ v.s. COZ v.s. KS v.s GS input nucleon DA 




0.0 0.2 0.4 0.6 0.8 



x B 



FIG. 7: Unpolarized cross section (in nb/sr) for backward 7*p — > pir° (upper panel) 
and for backward 7*p — > mr + (lower panel) as the function of Xb computed in the two 
component model for txN TDAs for Q 2 = 10 GeV 2 , u = —0.5 GeV 2 as a function of xb- CZ 
|| (solid lines), COZ (dotted lines), KS || (dashed lines) and GS @ (dash-dotted 
lines) nucleon DAs were used as inputs for our model. 

On the upper panel of Fig. we show the Q 2 dependence of the unpolarized differ- 
ential cross section of 7*p — > mr + for fixed £ = 0.25 which is characteristic for the J-lab 
kinematics and for A\ = 0. The plot exhibits the expected universal 1/Q 8 scaling behav- 
ior; the shape of tcN TDAs, indeed, affects only the overall normalization. On the lower 
panel of Fig. & similarly to |[49||, we show instead the same cross section as the function 
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of Q 2 for fixed W = 2.0 GeV and A^ = 0. Let us emphasize that the scaling behavior 
in the latter case is shadowed due to the fact that, for fixed W, the running of Q 2 also 
imposes variation of the scaling variable £. 



CZ v.s. COZ v.s. KS v.s GS input nucleon DA; £=0.25; A 2 T =0 




4 6 8 10 12 

Q 2 [GeV 2 ] 

CZ v.s. COZ v.s. KS v.s GS input nucleon DA; W=2.0 GeV; A 2 T =0 




0.001 [ 1 

4 6 8 10 12 

Q 2 [GeV 2 ] 



FIG. 8: Upper panel: unpolarized cross section (in nb/sr) for backward j*p — > nn + 
for fixed £ = 0.25 and = as a function of Q 2 in the two component model for ttN 
TDAs. Lower panel: unpolarized cross section (in nb/sr) for backward 7*p — > nir + 
for fixed W = 2.0 GeV and A^ = as a function of Q 2 in the two component model for 
ttN TDAs. CZ @ (solid lines), COZ g| (dotted lines), KS @ (dashed lines) and GS 
|43[| (dash-dotted lines) nucleon DAs are used as inputs for the model. 

Asymmetries, being ratios of the cross sections, are less sensitive to perturbative cor- 
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rections. Therefore, they are usually considered to be more reliable to test the factorized 
description of hard reactions. For the backward pion electroproduction, an evident can- 
didate is the transverse target single spin asymmetry (STSA) |^7| defined as: 

t-2-K \ / r2ir 



A=j^ QT d(p\M$\ 2 - J W d<p\M$\ 2> j (J ^ d(p\M% \ 2 ^j 
4^Im(X'(X)*) 



7T I7-I2 _ T IT/ 1 2 
l x I A/ 2 l X I 



(64) 



As argued in Sec. f|, within the two component model for ttN TDAs, the non-vanishing of 
the numerator in the last equality of (63) is achieved due to the interference of the spectral 
part contribution into ImZ(£) and of the nucleon pole part contribution into ReZ'(£). 

On Fig. we show the result of our calculation of the STSA for backward ir° and tt + 
electroproduction off protons for Q 2 = 10 GeV 2 and u = -0.5 GeV 2 . CZ §, COZ |3J , 



KS j|2] and GS |^3|| nucleon DAs are used as phenomenological input for our model. We 
conclude that STSA turns out to be sizable in the valence region. Its measurement should 
therefore be considered as a crucial test of the applicability of our collinear factorized 
scheme for backward pion electroproduction. 
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FIG. 9: Transverse target single spin asymmetry ( j64[) for backward 7*p^ — > pir (upper 
panel) and for backward 7*p"r — > nii + (lower panel) as a function of Xb computed with 
the two component model for tiN TDAs for Q 2 = 10 GeV 2 , u = —0.5 GeV 2 as a function 
of Xb- We show the results of our model with different input nucleon DAs: CZ (solid 
line), COZ (dotted line), KS (dashed line) and GS (dot dashed line) used as input. 



6. CONCLUSIONS 



For the first time, we have managed to build a consistent model of irN TDAs in their 
whole domain of definition. It satisfies general constraints imposed by the underlying 
QCD such as isospin symmetry, the Lorentz invariance manifested through the polynomi- 
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ality property of the Mellin moments of ttN TDAs in the light-cone momentum fractions, 
as well as the chiral properties. We used this model in the estimates of the unpolarized 
cross section and the transverse target single spin asymmetry for backward ir + and ir° 
electroproduction off protons. Our results make us hope for bright experimental prospects 
for measuring baryon to meson TDAs with high luminosity lepton beams such as COM- 
PASS, J-lab@ 12 GeV and EIC [fig] . Experimental data from J-lab@ 6 GeV on backward 
7r + , 7r°, T) and lo meson production are currently being analyzed [50]. We eagerly wait for 
the first evidences of the factorized picture for backward electroproduction reactions as 
suggested in our approach. 
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A. PARAMETRIZATION OF LEADING TWIST vriV TDAS 



The parametrization of the leading twist-3 ttN TDAs of given flavor contents suggested 
in ||T2[ which we employ in this paper reads: 

3 



A{P-nf 



-j-j- d\j 



3=1 



2tt 



e * Efc=i %k*k(P-n) l 



7r(p 7r )|O pTX (Ain, X 2 n, A 3 n)|A(pi)) 



= 6( Xl + x 2 + x 3 - 20 x i-^L [V? N {x u x 2 , x 3 ,£, A 2 )(PC) pT (PU) x 

+Al N (x h x 2 , x 3 , £, A 2 )(P 7 5 C) pr ( 7 5 Pf/) x + T? N (x u x 2 , x 3 , e, A 2 )(a Ptl C) pr {^PU) X 
+V? N (x 1 , x 2 , x 3 , f , A 2 )(PC) pr (AU) x + A* N { Xl ,x 2 , x 3 , £, A 2 )(P^C) pr ^ AU) X 



+T; N (x u x 2 , x 3 , f , A 2 )(a Pll C) pT (YAU) x + — Tf v (x u x 2 , x 3 , £, A 2 )(a PA C) pT (PU) 



ttN , 



AT 



+—T£ N (x 1 ,x 2 , x 3 , £, A 2 )(a PA C) pT (AU) J . 



(Al) 
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where is the pion weak decay constant and Jn is a constant, which determines the 
value of the dimensional nucleon wave function at the origin; U is the usual Dirac spinor 
and C is the charge conjugation matrix. We employ Dirac's "hat" notation: a = 7^a M and 
adopt the conventions: o^ v = |[7 M , 7 1/ ]; cr vu = v^a^, where is an arbitrary 4- vector. 



The relation of the parametrization (|A1| ) for nN TDAs to that of is given by: 



{Vx, A 1 ,T 1 y N \^ = (r^iVi, A,, T\} nN - j^{V 2 , a 2 , T 2 y N 

OW^Ig = {{V2,A 2 y N + l -{v x ,A x y N ) 



Jll & this paper ' 



& this paper 



rpirN I rpirN I I rpirN I 

' ; ' I 1 < ! I J 2 I Q & this paper ^ 2 1 I @] & this paper 

T'7r./V I / ^ rjiTvN 1 rrmN 1 rjinN cts-rriirN 

2 Ira = 1 2 Ti +T 2 +T3 _2 ^ T4 



rr.VI _ / 1 + ^rrmN , /-i , y-^nN 



T rlra = ( -^T + a + 0?? 
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Jl^ l & this paper 

(A2) 



|12l & this paper 



B. AN ALTERNATIVE FORM OF SPECTRAL REPRESENTATION FOR 
GPDS AND BARYON TO MESON TDAS 

1. GPD in the ERBL and DGLAP regions 

From (|lID, one can derive the following expressions for GPD in the DGLAP and the 
ERBL regions: 

1. for -1 < x < -f (DGLAP 1 region): 

2. For -f < x < £ (ERBL region): 

HM = j±- ^ dKF L,«^l^Y (B4) 

3. For f < x < 1 (DGLAP 2 region): 

= T 1 1 f 1 dtFU * 1 *®- 2 * ). (B5) 
s 1+1 s 
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2. Set of working formulas for irN TDAs in the ERBL-like and DGLAP-like 
regions 



In order to be able to compute irN TDAs from the spectral representation ([19]), we 
perform integrals over /i and A with the help of two 5-functions. We omit the index % 
referring to the choice of quark-diquark coordinates in the formulas of this Appendix. 

The resulting domain of integration in (k, 9) is defined by the inequalities: 

1 — k „ 1 — K 
-K/c<l; <9< ; 

-l+£ + 2w l-£ + 2w 

< K < 



—(w -2v + ] —^) < 6 < -- + — *— (™ + 2v + ] —^). (B6) 

2 1 + C 2 ' ~ ~ 2 1 + C 2 ' K 1 

Below, we summarize the explicit expressions for 7rN TDAs from the spectral repre- 
sentation flT9p in the ERBL-like and DGLAP-like regions. Let us introduce the following 
notation for the integrand: 

F (....) s F ( K , <L±^L, ^f^) ■ (B7) 

1. For w G [-1; -£] and u G [£'; 1 - £' + £] (DGLAP-like type I domain): 

1 — 2v4-w 1 — k. 

i r^+r- 



H(w,v,0 = —— dK d9F{....) (B8) 

2. For ioe[-l; -f] and v G [-£'; £'] (DGLAP-like type II domain): 

H{w,v,0 = Ti 7\i / rfK / d9F(....) (B9) 



3. For w G [-1; -£] and »e[-l + f- f ; -f] (DGLAP-like type I domain): 

H(w,v,£) = - — 2 I dK I d6F{....). (BIO) 



-i 



1 — K 



4. For w G [-^; f] and v G [£'; 1 - f + £'] (DGLAP-like type II domain): 

1 /-^rr" /-V 
H(w,v,0 = j^—^ dK d9F(....). (Bli; 
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5. For w G [-£; f] and v G [-£'; +£'] (ERBL-like domain): 

i / -f+ T i f (^+2«+i=i) 

H{w,v,t) = —— dK ddF(....). 

(l - o 2 i-±i±- y f -^^i^ 

6. For w G [-£ ; f] and v G [-1 + £ - -f] (DGLAP-like type II domain): 



;i-o 



1^— i: ~-*+^(t U +2t;+i=£) 



iZYtt;,i;,f) = - ^/ <*« / dOFL...). (B12) 

' 1 -l+g+2^ / 1-n 



i+e 



7. For w G [£; 1] and v G [— 1 — £+£'], the result coincides with ( Bll|) as it certainly 
should be, since this is the part of the same DGLAP-like type II domain. 

8. For w G [f; 1] and v G [£'; -£'] (DGLAP-like type II domain): 

1 ^ 



H{w,v^) = 7 ——r= dn d9F{....). (B13) 

(1 - £) 2 J-l+i+2v 



1+4 " 2 



9. For w 6 [£; 1] and f G [—1 +£ — £'; £'], the result coincides with ( |B12| ), since this 
is the part of the same DGLAP-like type II domain. 



C. ON THE RELEVANT GENERALIZED FUNCTIONS 

Sohotsky's formula (see e.g. Chapter II of |4(J) reads: 

= Tin5(x) + V- , (CI) 
x ± zl) a; 

where P stands for the Cauchy principal value prescription. The generalized function V-^ 
is then defined as -^V^ = —V\. For an arbitrary test function ip(x), 

P ±,, W W.^M^>. (C2) 



Employing flCl ) and ( p2[) one can establish the familiar relation: 



/I 1 

Ti^S'(x)-V^. (C3) 



dx x ±i0 x 2 

The formula ( |U3| ) concerns the conventional generalized functions dealing with the 
class of test functions defined on (— oo; oo) and sufficiently fast decreasing at the infinity. 
In our case, we have to consider a different class of generalized functions dealing with the 
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test functions defined on the interval [A; B] (A < 0, B > so that the singularity point 
x = belongs to the interval). Let <p(x) be a test function defined on the interval [A; B]. 
Then 



U '±(vl)M*)) = (4-(rhM*: 



dx \ x 



dx V x 



[A,B] 



x=B 

— lim 

x=A \JA 



+ 



dx 



X 



Let us consider the integral term at the r.h.s of eq. (C3): 



lim 

£^0 



+ 



B 



dx 



X 



X 



i<p(x) - <p(0)) 



x=B 



lim 

e-5-0 
B 



+ 



B 



dx 



^(y(x)-y(O)) 
x 



+ V dx 

x=A J A X 



1 <p{x) - cp(0) 



X 



Thus we conclude that 



(C4) 



(C5) 



x- 



X 



[A,B] 
x=B 



— (v- 



x=A J A x x 



[A,B] 



(C6) 



So finally we establish the following relation: 
1 



yxiCM) / [A,B] V x / [A,_B] 

w(o) + + p / ^- ?1 V U . 



(C7) 



that represents the version of (|C3|) adopted for the use on the finite interval. 



D. CALCULATION OF THE CONVOLUTION INTEGRALS 

In the calculation of Re/j ± ' ± ' ) (£) and Reljj' (£)> we encountered the following double 
principal value integrals 

l fc Ra p / r „ ! * l (^ i,( "' ,, ' ?,i (D1) 

"/! iw j^w T CZ t dv ¥hj) H ^ «■ (D2) 
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We propose here a strategy of computation of these integrals once H(w, v,£) is param- 



integrals when computing the real part of the elementary DVCS amplitude with GPDs 
parameterized through Radyushkin's factorized Ansatz. The following steps are to be 
performed: 

1. By interchanging the order of integration in ( pl|) and (|D2| ) , w and v integrals may 
be computed using the two delta functions. One is left with four integrations over 
the spectral parameters. 

2. After suitable change of variables, two principal value integrations can be performed 
analytically. 

3. The double integration over the remaining two spectral parameters is performed 
numerically. The corresponding integrands possess only logarithmic singularities 
which are perfectly integrable. In this way, we managed to reduce the problem of 
performing highly singular principal value integrals ( pi| ), ( |D2| ) to much less singular 
integration. This allows us to construct a stable and reliable numerical procedure. 

We present below the results for the principal value integral ( pi| ) for the case of the 
factorized Ansatz ( pOP with the profile h(fi, A), given by eq. (p4|). 



eterized with the help of the spectral representation ( |T9| ) with the use of the factorized 
Ansatz (PU[). The procedure generalizes the way of proceeding with the principal value 




l 



/_ 



dO Z 1 {a^ ± \n, 6, 0, & (±,±) (^ 9, 0, c (±) (^ 0) V(k, 9). 



(D3) 



Here rj is the sign factor: 



V 



(±.+) 



1: 



V 



(±-) 



-1 . 



(D4) 



The coefficient functions a^'^, b^'^ are defined as follows: 




(D5) 
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2 V 2 



&<+'+>(«, e, = ^-\(^-e) (i + O; 

&(+•-)(«, i) = 2e - i + a + o ■ (D6) 

The first (second) sign in the indices of t/ 1 * 1 '^, a ( - ± ' ± - ) , fe^'^ corresponds to that in the 
w ± £ (f ± £') denominator in ( [DT| ) respectively. The coefficient functions are defined 

as 

C W(k, = ~ (k(1 + + (1 - ± 2C) • (D7) 

The sign in the index corresponds to the one in the w ± £ denominator of 
The explicit expression for Zi(a,b,c) in (|D3|) reads: 



_ , , . a 3 , 9 3ca 2 6 2 a 3c 2 a 6 3 c 3 , 9 6 2 c 

_i(a, 6, c) = h 6a H 1 h be 

6 2 2 2 3 6 2 

a& 2 ac 2 \ / (6 — c ) 2 \ /a 2 c & 2 c\ / (a — 6) 2 



+ hr-^r log _ _____ u 



2 2 / ~° V a 2 J V 2 2 J to V c 2 

2\ /(„ h\2~ 



1 .1 /a 2 \ /(a -6) 
+a6c - log — log ' 



2 ° U 2 / ° \ c 2 
>og ( 1 - I) log - 2L l2 (5) + log ( I log ( - j + 2L, 2 ^1 - - 



where lAz^z) is the usual dilogarithm function 



(D8) 



Li 2 (z) = - I dz^-^ . (D9) 



One may check, that for real a, b and c, no imaginary part appears in ( P8|) as it should be, 
since it is the result of integration of a real function over a real interval. The imaginary 
part occurring from dilogarithms for | > 1 and | < is exactly canceled by the imaginary 
parts stemming from logarithms in the last line of ( |D8| ) due to the well known property 
of dilogarithm: 

Im(_i 2 (z + i0)) = -Trlog^ for z > 1 . (D10) 
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We now turn to the second principal value integral ( |D2| ) . For the factorized Ansatz 
(p0|), with the profile h(fj,, A), given by eq. (^4|) the result reads: 

i /-i-ie-e'i 



r i i 

V dw l 7V2 V d Vr ——H(w,v,0 

w,,/1 ' f dK I 2 de z 2 ( a { -' ±) (K, e, o, ft^fo 0: c(_) (^ 0) 

J Kn J — — =r^ 



(i-0 5 

60?/-^ 

(W) 



1-K 



+7T^ / rfK / / dd ^a^\n, 0, 0, ^(ic, £), c H («, 0) V(«, 0), 



-i 



(Dll) 

where t/ _,± ), a^ _,± \ &( _,± ), a* - * 1 are given by ( |D4|) - ( |D7|) and kq is defined by the equation 

C (-W) = 0: ^0 = ^^. (D12) 
The explicit expressions for Z 2 (a, b, c) and Z 2 (a, b, c) read 
Z 2 (a, b, c) 

1/ r Wr or \ / 7 \ n / / 7\2\ GIC (<2 lOg (<2 2 ) — Mog (& 2 ) ) 

= — (a - 6 + c)(5a + 36 + c) + (a - 6)clog ((a - by) - taV ' ta v " + 

2 a — o 

\ K ~ b 2 ) log {^P) + «6 log + 6c log (c 2 ) + a(6 - c) log 
--a{b + c) log — log' 



2 v y V & 2 / V c 2 



+ c, (, g (! - 5) l0g (Q - log (£) log ((! - £)') + 2L, (£) - 2Ll2 (!-£)); 
Z 2 (a, b, c) = Z 2 (a, b, c) + (^-^ + c (V - a 2 + afolog J . (D13) 

To complete the calculation of the real and imaginary parts of l\j (0), we also 
present the explicit expression for ( ^^f^ ~ ) anc l ( dJ< Iw"'^ I ( see ec L - (PI))- 

V / w=£ V / w=£ 



Using the formulas summarized in the Appendix |B 2j one may check that 



dH(w, t?,£) 
dw 

i f 1 ^ /" v Jfl m (-.,(0,1) Mi + - 2£ 0(1 + 



dn J 2 d9AV(K,9) j^ ' 1 ) 



i-e ' i-e 



■^ ( '^+o-^ga±|) , )i (D14) 



36 



where h^ 1,0 \^i, A) = j^h(fi, A); /i( 0,1 )(/i, A) = -§^h([i, A) and we use the fact that h van- 
ishes at the border of its domain of definition: h(—l,\) = 0. 
Finally, 

V dw J w=( 

= f dK J 1 d * ^)^3(a ( -' ±} («, 0, 0. 6 (_,:t) («» 0, ^(k, 0). 

(D15) 



where 

?! 
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Z 3 (a, 6, c) = (a - 6)(a + b + 2c) - a(6 + c) log ( ^ ) . (D16) 
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